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ABSTRACT 

We report the discovery of a very large, spatially extended Lya-emitting neb¬ 
ula at z = 2.656 associated with a luminous mid-infrared source. The bright 
mid-infrared source (F 24 /rm = 0.86 mJy) was first detected in observations made 
using the Spitzer Space Telescope. Existing broad-band imaging data from the 
NO AO Deep Wide-Field Survey revealed the mid-infrared source to be associ¬ 
ated with a diffuse, spatially extended, optical counterpart in the band. 
Spectroscopy and further imaging of this target reveals that the optical source is 
an almost purely line-emitting nebula with little, if any, detectable diffuse con¬ 
tinuum emission. The Lya nebula has a luminosity of Li^ya ~ 1.7 X 10^^ erg s ^ 
and an extent of at least 20" (160 kpc). Its central ~ 8" shows an ordered, 
monotonic velocity prohle; interpreted as rotation, this region encloses a mass 
M K, Qx IO^^Mq. Several sources lie within the nebula. The central region of the 
nebula shows narrow (~ 365 km s“^) emission lines of C IV and He II. The mid- 
infrared source is a compact object lying within the nebula, but offset from the 


^National Optical Astronomy Observatory, 950 N. Cherry Ave., Tucson, AZ 85719; dey@noao.edu 
^California Institute of Technology, Pasadena, CA 91125 

^Spitzer Science Center, California Institute of Technology, MS 220-6, Pasadena, CA 91125 
'^Dept. of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544-1001 
®W. M. Keck Observatory, 65-1120 Mamalahoa Hwy., Kamuela, HI 96743 
® Steward Observatory, University of Arizona, Tucson, AZ 85721 
^Astronomy Department, Cornell University, Ithaca, NY 14853 

®Jet Propulsion Laboratory, California Institute of Technology, MC 169-327, 4800 Oak Grove Drive, 
Pasadena, CA 91109 



2 


center by a projected distance of ~ 2'/5 (20 kpc), and likely to be an enshrouded 
AGN. A young star-forming galaxy lies near the northern end of the nebula. We 
suggest that the nebula is a site of recent multiple galaxy and AGN formation, 
with the spatial distribution of galaxies within the nebula perhaps tracking the 
formation history of the system. 

Subject headings: galaxies:formation-galaxies:high-redshift-galaxies:starburst 


1. Introduction 

How do the most massive galaxies form? This is a question whose answer may provide 
one of the central tests of hierarchical galaxy formation theories. Our current theoretical 
paradigm suggests that the most massive galaxies form, like the rest of the galaxy population, 
by the hierarchical merging of smaller units and that this process occurs through the history 
of the universe. In seeming contrast, however, three lines of evidence suggest that at least 
some of the most massive galaxies may have had more dramatic origins. First, observations 
of the most massive galaxies at redshifts out to 2 : < 1.5 suggest that the bulk of the stars in 
these objects formed at much higher redshifts in a brief time period characterized by very 
high star-formation rates (> 300Mq yr“^) (e.g., McGarthy et ah 2004). Second, the host 
galaxies of powerful radio sources, which appear to be dynamically relaxed massive ellipticals 
at redshifts out to one (e.g., Zirm et al. 2003), at higher redshifts ( 2 :^ 3 ) are surrounded 
by large (100—200 kpc), luminous (> lO^^ergs"^) Lya nebulae and small, embedded, star¬ 
forming objects (e.g., Reuland et al. 2003). This suggests that these host galaxies are formed 
in a rather spectacular gravitational collapse event rather than by the slow build-up of many 
merging events over time (Reuland et al. 2003). The discovery of large, diffuse Lya haloes in 
several surveys and their association with large galaxy overdensities adds further support to 
the idea that these haloes mark the formation sites of massive galaxies (e.g.. Keel et al. 1999; 
Steidel et al. 2000; Francis et al. 2001; Palunas et al. 2004; Matsuda et al. 2004). Third, 
the recent discoveries that dust-enshrouded galaxies selected in sub-mm surveys have very 
large star-formation rates, large correlation length and high median redshift suggest that 
this population may represent the formation epoch of massive galaxies (e.g., Blain et al. 
2004; Ghapman et al. 2004). Taken together, the empirical picture that emerges for massive 
galaxy formation is one of a violent, perhaps brief epoch of dust-enshrouded star-formation 
followed by a more gradual unveiling of the galaxy by the dissipation of its dusty, gaseous 
envelope. 

It is clearly of great interest, therefore, to hud and characterize the formation sites of 
massive galaxies. Thus far we have been limited to studying the (>100 kpc) Lya halos 
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and environments of luminous radio galaxies, which may be affected in peculiar ways by 
the presence of the active galactic nucleus and its ejecta, and to the similarly large Lya 
nebulae in two regions: four associated with a massive structure at z = 2.38 (Francis et ah 
2001; Palunas et ah 2004) and two others associated with a structure at z = 3.09 (Steidel 
et al. 2000; Matsuda et ah 2004). It is worth noting that the only narrow-band emission 
line surveys that have yielded discoveries of these spectacular Lya haloes have been targeted 
at large mass overdensities, so their mutual association may be misleading. Ground-based 
sub-mm surveys provide a less biased approach in principle, but so far have been limited 
by technology: the current suite of bolometer arrays provide poor positional accuracy and 
limited helds of view which are not suited to wide-held surveys. 

In this paper, we present observations that suggest that mid-IR surveys with Spitzer 
Space Telescope may provide an alternate approach: we report on observations of a luminous, 
very spatially extended Lya nebula at z = 2.6562 discovered by virtue of its associated 24/rm 
emission. We present the description of our observations of this nebula in §2 and our results 
and measurements in §3. In §4, we discuss the implications of this discovery on future studies 
of massive galaxy formation. 

Throughout this paper we use Hq = 70 km s“^ Mpc“\ Gm = 0.3, Ga = 0.7. At a red- 
shift = 2.6562 the luminosity distance is di = 21.946 Gpc and 1 " corresponds to a physical 
size of 7.96 kpc. 


2. Observations 

The optical photometry of this source is derived primarily from the images obtained by 
the NOAO Deep Wide-Field Survey (NDWFS^; Jannuzi & Dey 1999; Jannuzi et al. 2005; 
Dey et al. 2005). The NDWFS is a deep, ground-based, optical and near-infrared imaging 
survey of two 9.3 square degree helds, one in Bootes and one in Getus, conducted using 
the 4m and 2.1m telescopes of the National Optical Astronomy Observatory. The survey 
is designed to study the evolution and clustering of galaxies over a wide range in redshift, 
and reaches median 5a point-source depths in the Bw,R,I and K bands of ~27.1, 26.1, 
25.4 and 19.0 mag (Vega) respectively. The data products for the NDWFS will be described 
elsewhere (Dey et al. 2005; Jannuzi et al. 2005, in preparation). The specihc optical image 
stacks used in this paper have R, and / band exposure times of 2.67, 2.67 and 3.33 
hours and reach 5a point-source depths of 27.31, 26.23 and 25.47 mag respectively. Stars in 
all three image stacks have FWHM^l". The source is undetected in the near-infrared Ks 


^See also http://www.noao.edu/noaodeep/ 
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data of the NDWFS. The NDWFS astrometry is tied to the reference frame dehned by stars 
from the USNO A-2 catalog (Monet et ah 1998). 

The NDWFS Bootes held was mapped at 24, 70 and 160pm using the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et ah 2004) in late January 2004 to la rms depths 
of SlpJy, 5mJy and 18mJy respectively. The data were reduced by the MIPS GTO team. 
Details of the survey, including the mapping strategy, data reduction and the resulting 
catalog, will be discussed elsewhere (LeFloc’h et al. 2005, in preparation). The entire 
Bootes held was also mapped at 3.6, 4.5, 5.6 and 8.0pm using the Infrared Array Camera 
(IRAC; Fazio et al. 2004) on Spitzer, details of the IRAC observations may be found in 
Eisenhardt et al. (2004). Both the MIPS and IRAC data are astrometrically calibrated to 
the 2MASS astrometric frame. 

Initial optical identihcations of the MIPS sources with 24pm hux densities brighter than 
0.75mJy were done using the image catalogs from the NDWFS. The MIPS 24pm sources are 
all unresolved at the 6" angular resolution of the Spitzer Spaee Telescope. The images of all 
MIPS sources with optical associations fainter than 1=23 mag were visually inspected. A 
24pm source SST24 J1434110+331733 with a hux density of 0.86mJy was found to lie within 
a faint, dihuse, spatially extended optical source with blue optical colors, which suggested 
that the colors may be the result of strong line emission in the Bw band. 

Spectroscopic observations of the MIPS 24pm candidate SST24 J1434110+331733 were 
obtained using the Low Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) on the 
Keck I telescope on U.T. 2004 May 21 and 2004 June 22. The hrst set of observations, 
four 30 min exposures, were obtained in 0.8" seeing through a 1.0" wide and 11" long slitlet 
oriented at PA=0°. LRIS was conhgured with the 5600A dichroic, the 4001/mm (Ac=3400A, 
l.OOA pix“^) grism on the blue side and the 4001/mm (Af,=8500A, 1.86A pix“^) grating on 
the red side. A second set of hve 30 min observations were obtained in 1" seeing through 
a 1.5" wide and 30" long slitlet oriented at PA=19.4°, with the 6800A dichroic, 3001/mm 
(Ac=5000A, I. 43 A pix“^) grism, and the same red grating. In each observation set, the tele¬ 
scope was offset by a few arcsec between individual 30 min exposures. Relative spectropho- 
tometric calibrations were performed using observations of Wolf 1346 (Oke 1990; Massey 
et al. 1988; Massey & Cronwall 1990). Figure 1 shows the location and orientation of the 
LRIS slitlets on the Bw image of the source. 

Additional optical observations of the held were obtained using the MiniMOSAIC cam¬ 
era (“MiniMo”) on the WIYN 3.5-m telescope of the Kitt Peak National Observatory on the 
nights of U.T. 2004 June 14 ~ 17. A total of six 1200 s images were obtained through the 62A- 
wide narrow-band Liter centered at 4448A (hereinafter NB4448, also called ‘W20” in the 
WIYN Liter set; centered on the redshifted Lya line) and three 900 s images were obtained 
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2004 May 21 PA=0 



2004 June 22 PA=19.4 



Fig. 1.— The left and right panels show the location, orientation and size of the LRIS 
slitlets superimposed on the Bw image of SST24 J1434110+331733. The origin corresponds 
to 14^34™l(h977 + 33°17'3(y.'87, the centroid of the Lya emission as determined from the 
narrow-band image (please see text). Note that only the slitlet at PA=19.4 (on the right) 
includes light from galaxy B. 
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through a Gunn G hlter. The seeing varied from night to night, but point sources in the 
hnal NB4448 and G-band image stacks have an effective FWHM of &!7 and 0'/8 respectively. 
The data were calibrated using observations of the spectrophotometric standard Feige 110 
(Oke 1990; Stone 1977) through the NB4448 filter and Landolt field SA112 (Landolt 1992) 
through the Gunn G filter on photometric nights. 

All the WIYN optical imaging data were astrometrically calibrated using the USNO 
A-2 catalog and tied to the same astrometric frame as the NDWFS imaging data. There is a 
small astrometric offset between the IRAG and MIPS 2MASS frame and the NDWFS/USNO 
A-2 frame of roughly Cf.'S which we corrected before comparing the Spitzer and ground-based 
data. All positions quoted here are on the NDWFS/USNO A-2 reference frame. 


3. Results 

3.1. Morphology and Astrometry 

Image cutouts of the B]y, G, R and / images are shown in figure 2, with the contour 
map from the NB4448 image overlaid on the /-band image^. The system is diffuse and 
relatively bright in the Bw band, which is largely dominated by the Lya line emission; there 
is a fairly compact galaxy at its northern extent which is labelled ‘A’ in figure 2. In contrast 
to the Bw image, the G, R and / band cutouts show only little (if any) diffuse emission, bnt 
do show some light associated with the regions ‘A’ and ‘B’. The Lya emission appears to 
avoid ‘A’. The center of the Lya emission lies near ‘B’, but slightly to its east; the spectrum 
of this central region is therefore contaminated by continuum emission from ‘B’. 

The NDWFS astrometry is tied to a frame defined by moderately bright (12—17 mag) 
stars in the USNO A-2 catalog, whereas the MIPS astrometry is calibrated using 2MASS. 
Roughly 98% of the bright (F 24 /im > 0.75 mJy) MIPS sources have well determined optical 
connterparts in the NDWFS data, and comparisons of the relative astrometry between these 
MIPS and NDWFS matches show a small, but significant relative offset of Aa = 0738 and 
AS = 0714 between the two frames. The offset arises because of a small magnitude-dependent 
astrometric error internal to the USNO A-2 catalog (Monet et al. 2003). The MIPS source is 
located at Id^dd^lO^OSS -|- 33°17'32757 in the J2000 2MASS reference frame. In the optical 
reference frame of the NDWFS, the source position is Id^dd^'lO'lOSS -|-33°17'32771 (J2000). 
Taking this offset into account, the MIPS source is found to he on the northern end of the 


^Larger finding charts of the field and astrometry of offset stars may be determined from image cutouts 
of the NDWFS available from the NOAO Science Archive (http://www.archive.noao.edu/ndwfs/). 
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Fig. 2.— Image cutouts of the SST24 J1434110+331733 in G-, R- and /-bands. The 
origin (0,0) in all panels is at Id^Od'^lFOOS -|- 33°17'29'.'84 (J2000). The red crosses in the 
four panels mark the location of the 24/rm source {Bw panel), the 3.6/rm source {G panel), 
the IRAC 4.5/im source {R panel) and the 8.0/rm source (/ panel). The position of the 
source in the MIPS 24/rm and IRAC bands agrees to < 075. The contours in the I panel 
denote the emission detected in the NB4448 hlter centered on the Lya emission line. The 
continuum emission in this hlter has not been subtracted, but contributes less than 10% of 
the total hux at all radii. The contours drawn correspond to line hux surface brightnesses 
of approximately [0.2,0.4,0.6,0.8] xl0“^®erg s“^ cm“^ arcsec”^. 





















diffuse emission observed in the NDWFS Bw and NB4448 images, just to the southwest of 
galaxy A (£g 2). The location of the MIPS source is marked by a cross, and is believed to 
be accurate to ~ 0'/5. It lies ~ 0'/8 west and T/6 south of galaxy A and ~ T/1 east and 2'/6 
north of galaxy B. The Lya emission is centered roughly 2'.'5, or 20 kpc (projected), away 
from the MIPS source. 

The position of the IRAC source detected in the 3.6/im, 4.5/im, and 8/im bands agree 
with the MIPS 24/im source position to <0'.'5, and confirm the relative (non-central) location 
of the mid-infrared source within the Lya nebula. 


3.2. Photometry 

Table 1 presents our photometric data on SST24 J1434110-I-331733. The optical pho¬ 
tometry was measured in a large aperture of radius 10" centered on the centroid of the 
Lyo halo. The flux in the By/ band is dominated by the Lya emission, and the optical 
measurement aperture for all bands includes galaxies A and B. These optical measurements 
should therefore be strictly treated as upper limits to the continuum emission from the re¬ 
gion. Higher spatial resolution imaging and spectroscopy along PAr^OO® will be necessary 
to disentangle the contribution of galaxy B and determine the true continuum contribution 
from sources within the nebula. 


3.3. Redshifts of the Nebula and Coincident Sources 

The spectroscopic observations reveal a very complex system within the extent of the 
nebula, comprised of several continuum sources. Figure 3 shows the two-dimensional spec¬ 
trum recorded on the blue side of LRIS in PA=19.4°, and figure 4 shows the extracted 
one-dimensional spectrum of the system obtained through a (wide) 10'.'9xT.'5 aperture from 
these data. Emission line measurements are presented in table 2. The redshift of the system 
is well established ai z = 2.6562 through the detection of Lyo;, C IV, He II and weak C III]. 

There are three main continuum contributors to this large-aperture spectrum: (i) galaxy 
A (figure 5), (ii) a source T.'6 southwest of galaxy A along the PA=19.4° slitlet which may 
be associated with the MIPS source (figure 6), and (iii) a source 477 southwest of A, and 
close to the location of galaxy B (figure 7). We discuss each of these in turn. 

The spectrum of galaxy A (figure 5) reveals a “Lyman-break galaxy” at a redshift 
of X ~ 2.656, characterized by a flat F7 continuum, interstellar absorption lines and no 
detectable Lya emission. The interstellar absorption lines observed in typical Lyman-break 
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Fig. 3.— Two-dimensional spectrum of the Lya nebula associated with SST24 
J1434110-I-331733 obtained through a I'.'S wide slitlet oriented in PA=19.4° and using LRIS- 
B on the Keck I telescope. The three continuum streaks and main emission lines are labelled. 
Please see text for details. 


Table 1. Photometry of SST24 J1434110-I-331733 


Band 

Wavelength 

(^tm) 

Width 

(/im) 

Aperture Center 
(J2000) 

Aperture 

Radius 

Magnitude 

(Vega) 

F. 

(erg s~^ cm~^ Hz~^) 

NDWFS Bw 

0.422 

0.112 

14'‘34’"ia;977 +33°17'30'.'87 

10" 

22.64 ± 0.05 

3.18 ±0.15 X 10-29 

Gunn G 

0.492 

0.067 

” 

10" 

23.02 ± 0.20 

2.39 ±0.45 X 10-29 

NDWFS R 

0.659 

0.183 

” 

10" 

22.24 ±0.08 

3.92 ±0.30 X 10-29 

NDWFS I 

0.808 

0.172 

” 

10" 

22.02 ± 0.15 

3.80 ±0.53 X 10-29 

NDWFS Ks 

2.167 

0.295 


6" 

> 18.12=" 

< 3.67 X 10-2®=" 

IRAC 3.6/^m 

3.56 

0.75 

14'*34’"10?979 + 33°17'32'.'48 

6" 

17.83 ± 0.13 

2.04 ±0.24 X 10-2® 

IRAC 4.5/^m 

4.52 

1.01 

” 

6" 

17.09 ± 0.13 

2.63 ±0.31 X 10-2® 

IRAC 5.8/^m 

5.73 

1.42 

” 

6" 

15.93 ± 0.36 

4.97 ± 1.66 X 10-2® 

IRAC 8fim 

7.91 

2.93 

” 

6" 

14.79 ± 0.21 

7.69 ± 1.49 X 10-2® 

MIPS 24/im 

24.0 

5.3 

14'‘34’"10?988 +33°17'32771 

b 

9.83 ± 0.06 

8.56 ±0.51 X 10-27 

MIPS TO^m 

70.0 

19.0 




< 2.5 X 10-25 

MIPS 160/.tni 

160.0 

34.5 




< 9.0 X 10-25 

WSRT 1400MHz'= 

20cm 


14'‘34’"1K78 +33°17'2774 



1.4 ±0.4 X 10-27 


^The Ks upper limit is 2a in a 6" radius aperture. 

*^The MIPS photometry was determined using PSF-fitting techniques. 

"^WSRT measurement is from the deep 1400MHz survey of de Vries et al. (2002). The source is unresolved in the 13^^x27'' beam. 
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galaxies can be blue shifted due to winds, and the redshift estimate based on these is therefore 
only good to ~ 1000 km s“^. The absorption line spectrum seen in hgure 4 is almost entirely 
due to the continuum contribution from galaxy A to the large aperture extraction. 

The spectrum of the second continuum source (hgure 6) is quite red and shows weak 
C IV and C III] emission at z = 2.656, i.e., the redshift of the nebula. Its proximity to 
the MIPS 24/im location suggests that it may be a star-forming galaxy reddened by dust. 
There are some possible continuum features visible in the spectrum of this source, that are 
consistent with a young, reddened star-forming object. No obvious emission lines uniquely 
due to an AGN are detected in this source. 

The spectra of galaxy A and the MIPS source both show a fairly broad absorption 
feature at A5720A. The signal-to-noise ratio in this region is low, and sky subtraction using 
traditional htting techniques on very extended faint objects can sometimes result in such 
features. If this feature is real, it is most likely to be due to a foreground absorption system, 
perhaps arising from either MgIIAA2800 at z zz 1.04 or an Fell 2600A absorption complex 
at ~ 1.20.” 

The third continuum source (hgure 7) is the most intriguing: it lies near the center 
of the Lya nebula and apparently shows strong associated C IV and He II emission at the 
redshift of the nebula. This continuum source also shows an emission feature at 5081. OA, 
the identihcation of which is uncertain. This feature is unlikely to be an emission line at the 
redshift of the Lya nebula, since it does not match any strong features normally observed in 
high-redshift galaxies or AGN: at a rest wavelength of 1389. oA, it could only be associated 
with blue-shifted SiIVAA1393.8,1402.8 emission, and would be the only feature with such a 
blue shift relative to the He II redshift, which seems very improbable. This spectral line 
also does not appear in the PA=0° spectrum, strongly suggesting that it is associated with 
galaxy B. Since the 5081.6A line does not correspond to any normal species at the redshift 
of the Lya nebula, we therefore suspect that galaxy B is an interloping system that does not 
he at the redshift of the nebula. Given the asymmetry of the 5081.6A emission line and the 
break in the continuum across the line, we suggest that galaxy B is a background source at 
a redshift of 3.180. Galaxy B is located too far away from the centroid of the MIPS 24/im 
emission and is therefore also not the source of the infrared emission. The spatial coincidence 
of the 5081.6A feature and the He II and G IV emisson lines in our spectrum must therefore 
be due to the unfortunate circumstance of the two contributing sources, the emission line 
source at z = 2.656 and the background galaxy at z = 3.18, being located approximately 
along the short-axis of our slit, parallel to the dispersion direction. It is also likely that 
the continuum emission that we observe at this location is primarily due to galaxy B, and 
unrelated to the nebula. 
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The source near the top of the slit in hgure 3 shows an asymmetric emission line at 
5188.SA and a break across the line; this emission line is consistent with Lya at a redshift 
of 3.268, and this galaxy is also a backgronnd source. 

We are left, therefore, with the interpretation that two nnrelated backgronnd sonrces are 
present in this l'.'5x30" slitlet, both at redshifts above three. We are conhdent the emission 
line detections associated with these objects are not instrnmental effects, since they repeat 
on all the individual exposures and move in concert with their associated continnnm when 
the targets were offset along the slit. It is not exceptional to hnd a serendipitous source with 
z > 3 during LRIS spectroscopy; onr experience is that one snch sonrce appears in abont 
15% of the masks used. Given the snrface density of x ~ 3 Lyman break galaxies {zzl.2 per 
sq. arcmin for R < 25.5 mag; (Steidel et ah 2003)), the Poisson probability of hnding 2 in 
a single mask is 1.8%, and in a single slitlet is 0.01%. Note that these are lower limits to 
the probabilities since z ^ 3 Lyman break galaxies are clustered. The redshift separation 
between the two backgronnd galaxies snggests a minimnm comoving relative distance of 
~80 Mpc. Further spectroscopy of this region is clearly of interest in order to nnderstand 
the held galaxy redshift distribntion. 


3.4. Spatial Extent of the Lya Nebula 

A visnal examination of the NB4448 image at a hard stretch shows that the high snrface 
brightness portion of the Lyo; nebula is at least 10" in size, bnt that there is faint diffnse 
emission extending over larger scales. In order to explore the fnll spatial extent of the nebula, 
we compared its radial prohle with that of stars measured in the same image. Fignre 8 shows 
the integrated radial prohle of the nebula compared with the median prohle determined from 
10 stars in the NB4448 image, normalized to the same total hux. It is clear that the nebula 
is spatially extended to a radius of at least 15" (nearly 120 kpc)! There are hints that 
the nebula may be even more extended, but deeper narrow-band imaging observations are 
necessary to conhrm this. 

The emission seen in the narrow-band hlter is not signihcantly ahected by continnnm 
emission. The Gnnn G hlter only spans the wavelength range 1240—1460A in the rest frame 
of the system, and therefore excludes the redshifted Lya, G IV and He II line emission from 
the nebula. This band provides a good measnrement of the line-free hnx from the region 
which we can use to estimate (and subtract) the continuum contribution in the NB4448 hlter. 
In an apertnre of 10" radius, the continuum contribution is only 0.225 x 10“^® erg s“^ cm“^ in 
the NB4448 hlter, i.e., only 7% of the total hux. This estimate assumes that the continuum 
components are dne to interlopers and contribnte across the entire wavelength region of 
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Fig. 4.— One-dimensional spectrnm of the Lya nebnla associated with SST24 
J1434110-I-331733 smoothed by lOA. The apertnre is I'.'S wide and Kf-'Q long and oriented 
at PA=19.4°. The red line shows the noise spectrnm offset by — 1/iJy. The absorption lines 
arise predominantly from galaxy A, with some contribntion from galaxy B noticeable at 
A > 5100A. The interloping emission line at A=5081.6A is only observed in the PA=19.4° 
observation, and is therefore thonght to arise from galaxy B, which lies just to the west of 
the emission center of the C IV and He II line emission. 
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Fig. 5.— One-dimensional spectrum of Galaxy A extracted in an aperture I'.'S x 0"95 
from the PA=19.4° observation. The dashed line overlaid on the spectrum is the composite 
spectrum of a ^ ~ 3 star-forming galaxy constructed from a sample showing strong Lya 
absorption (Shapley et ah 2003). 
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Fig. 6.— One-dimensional spectrum of the continuum object closest to the location of the 
mid-infrared source. The extraction aperture is 1"5 x 0'.'95 from PA=19.4° observation. The 
spectrum is much redder than that of galaxy A (figure 5) but shows weak C IV. 
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Fig. 7.— One-dimensional spectrum of central region and galaxy B from the PA=19.4° 
observation. The extraction aperture is F.'5 x Cf.'QS. The spectrum cannot be disentangled 
cleanly since the two contributing components lie parallel to the dispersion direction, but 
we suggest that the 5081.bA line and the continuum emission are dominated by galaxy B, 
whereas the Lya, C IV and He II emission arise from the central region of the nebula. 
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the narrow-band filter. If we instead assume that the continuum arises from sources at the 
redshift of the nebula, then the continuum contribution within the narrow-band hlter would 
be reduced due to Lya forest absorption in the intergalactic medium. Using the Madau 
(1995) approximation for the total continuum absorption, we estimate that the continuum 
contribution would be about 13% smaller under this assumption, and the total derived Lya 
flux imperceptibly larger. 


3.5. Velocity Structure in the Nebula 

Our spectroscopic observations are not sufficiently deep to detect the nebula across its 
full extent, and we restrict our discussion here to the central 8" region. 

Lya is a resonance line, and its profile is generally dominated by radiative transfer 
effects, being affected by resonant scattering and absorption by gas and dust. It is therefore 
generally a poor tracer of the gravitational potential or the physical properties of the gas. 
However, at present, this is the only real tracer we have of this system; we therefore use it 
to provide some indication of the kinematic properties of the system with these caveats in 
mind. 

We extracted spectra in 0'.'4 (3-pixel) apertures from each of our spectroscopic observa¬ 
tions and £t the Lya profile with a Gaussian model. The fits are fairly good approximations 
to the emission line profile; the profiles are only very slightly asymmetric across most of the 
nebula. We also performed similar fits on the [OI]A5577 night sky emission line, and found 
that the tilt this line varied by < ±5 km s“^ across the portion of the slitlets covering the 
Lya nebula. 

Figure 9 shows the velocity structure in the Lya line as determined by these fits. The 
seeing for both observations was typically twice the aperture, and neighboring points in the 
figure are therefore correlated. The nebula exhibits a monotonic change across the central 
k.A" in both the central velocity and the velocity dispersion, with deviations from this trend 
only near the extrema where the signal-to-noise ratio is low. 

The remarkable aspect of the measured profile is that its velocity varies so uniformly 
across the central region. Fitting a straight line to the data from the two observations using 
a robust least absolute deviation algorithm results in a velocity profile across this central 
region of n = 36.2 — 124.Ir km s“^ in PA=0° and v = 9.7 —75.5r km s“^ in PA=19.4°, where 
r is the angular separation in arcseconds. This smooth variation is unlikely to be due to a 
smooth change in the absorption by neutral hydrogen, since the line profile does not show 
strong evidence for this. If it is due to the gas kinematics, then the variation could be due 
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Fig. 8.— The integrated flux of Lya as a function of radius as measured in the NB4448 
image (filled circles). The solid line shows the mean stellar profile in the image, determined 
using 10 stars. The filled circles represent the integrated Lya profile after subtraction of the 
continuum estimated using the Gunn G-band observations. The dashed vertical line at a 
radius of 10'' denotes our chosen (conservative) aperture size for most of the measurements 
presented in this paper. The horizontal dashed line is for reference. 
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to infall, outflow or rotation, and it is not easy to discriminate between these options with 
the present data. 


3.6. Spectral Energy Distribution 

Figure 10 shows the spectral energy distribution (hereinafter SED) of SST24 J1434110+331733. 
The optical data points should be regarded as upper limits, since they include continuum 
contributions from galaxies A and B. 

A model for the SED of Arp 220 (constructed from the models of Silva et ah (1998) and 
Spoon et ah (2004)) placed at z = 2.6562 and scaled by a factor of 26 is shown for comparison 
in figure 10. The scaling factor was chosen to match the observed IRAC 8/um flux density of 
SST24 J1434110+331733. The SED of Arp 220 differs from that of SST24 J1434110+331733 
in two ways. First, the rest-frame UV emission observed in SST24 J1434110-I-331733 is much 
larger than that from Arp 220. This is not surprising: as mentioned before, the UV points 
are contaminated by other sources (and the bluest point by the strong Lya line emission). 
Second, the 5.8/im and 8.0/im IRAC points lie signihcantly above the prediction of the 
Arp 220 template. This region of the spectrum is of interest because it is where the SED of 
Arp 220 tranisitions from being dominated by starlight at shorter wavelengths to being dom¬ 
inated by dust emission at longer wavelengths. The SED of SST24 J1434110-I-331733 shows 
no evidence for such a transition, suggesting that perhaps there is little or no signihcant 
contribution from starlight in this portion of the spectrum. Indeed, the mid-infrared con¬ 
tinuum measurements from Spitzer (3.6/im — 24/im) are well-modeled by a simple quadratic 
power-law, fiy oc z/“^, also suggesting that the Spitzer bands are not strongly affected by 
emission from polycyclic aromatic hydrocarbons (PAHs), of which the strongest features are 
redshifted out of the IRAC bandpasses. 

The mid-infrared color of SST24 J1434110-I-331733 lies in the region of parameter space 
more generally occupied by quasars than by normal star-forming galaxies (e.g.. Stern et al. 

2005; Lacy et al. 2004). Figure 10 also shows the SED of the Seyfert 1 galaxy Mrk 231 (e.g., 

Ivison et al. 2004) placed at the redshift of SST24 J1434110-I-331733 and scaled by a factor of 
2. The SED of Mrk 231 is a better fit to the observed 3.6/im to 24/im photometry of SST24 
J1434110-I-331733, and suggests that this source is an enshrouded AGN. The observed excess 
in the MIPS 24/im band over the Mrk 231 SED may be due to 7.7/im PAH emission. 

The most likely emission mechanism responsible for the mid-infrared luminosity is dust 
emission. Emission at these wavelengths would require dust at warm temperatures (> 500 
K), higher than those typically seen in normal or ultraluminous star-forming galaxies. How- 
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Fig. 9.— Kinematics of the gas in the Lya nebula surrounding SST24 J1434110+331733. 
The top two panels show the spatial variation in central velocity (left) and velocity dispersion 
(right) determined from Gaussian hts to the Lya emission line prohle measured using the 
PA=0° observation. The lower panels show the measurements from the PA=19.4° observa¬ 
tion. The origin for the ordinate is dehned as the spatial location of the centroid of the Lya 
emission (i.e., the same as in hgure 1). The velocity origin is dehned using the wavelength 
of the He IIA1640 emission line. 
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Fig. 10.— The spectral energy distribution of SST24 J1434110+331733. The optical points 
(Aobs < l/iui) should be regarded as upper limits, since they include continuum contributions 
from galaxies A and B. The red dashed curve is the SED of Arp 220 (A < 5/rm from Spoon 
et ah (2004) and A > S/im from Silva et ah (1998)), scaled in luminosity by a factor of 26, 
for comparison. Note the discrepancy between the two objects near Aobs ~ 10 —20/rm, where 
the SED of Arp 220 transitions from being dominated by starlight to being dominated by 
dust emission; SST24 J1434110+331733 shows no obvious evidence for such a transition. 
The blue solid line represents the SED of Mrk 231 from Ivison et ah (2004), scaled by a 
factor of 2. 
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ever, a power-law spectrum of oc can also be produced by synchrotron radiation from 
power-law distribution of electrons, with N{E) oc E~^. At longer wavelengths, synchrotron 
self-absorption sets in and the spectrum should take on a shape. We can use the faint ra¬ 
dio detection of 0.14 mJy and estimate the turnover frequency to be ~ 1.2 x 10^^ Hz (2.5mm). 

This is not a very useful constraint, especially since a synchrotron spectrum of this sort has 
never been observed in an extragalactic source. If the slope of the radio spectrum is similar 
to that of Mrk 231 and its extrapolation to the mid-infrared signihcantly undershoots the 
observed mid-infrared flux densities, then the emission mechanism for the mid-infrared is 
thermal emission from dust grains. Observations in the sub-mm to radio wavelengths will 
be necessary to better constrain the radiation mechanism for the IR portion of the SED. 

If SST24 J1434110-I-331733 has a similar SED to that of Arp 220 (Mrk 231), the compar¬ 
ison suggests that it has a total far-infrared luminosity of Lir < 7.5 X 10^3^0 (1.1 X 10 ^ 3 ^ 0 )! 

We emphasize, however, that these estimates are very uncertain at present, since the bulk of 
the luminosity emerges at wavelengths at which the flux density of SST24 J1434110-I-331733 
is currently unknown and require a large bolometric correction. Measurements at sub-mm 
wavelengths are clearly crucial to establishing the shape of the SED near its peak and pro¬ 
viding a better measure of the bolometric luminosity. If the SED is similar to that of Arp 220 
(Mrk 231), we would expect 450/im and 850/im flux densities of ^120 mJy (8 mJy) and 
<40 mJy (2 mJy) respectively. Clearly, these predictions are all within the reach of the 
current sub-mm observations. 

We note that the mid-infrared-to-radio flux density ratio observed in SST24 J1434110-I-331733, 
g 24 = log[/i,(24/im)//i,(20cm)] ~ 0.8, shows remarkable similarity to that of starburst galax¬ 
ies, which typically have g 24 = log[/i/(24/4m)//i.(20cm)] = 0.84±0.28 (Appleton et al. 2004). 

Mrk 231 also falls on this relation, although it is an obscured Seyfert 1; this is because its 
mid-infrared and radio emission are dominated by flux from the circumnuclear starburst 
rather than the active nucleus. 


4. Discussion 

In the previous section, we presented our various observations of SST24 J1434110±331733. 
Here, we discuss the implications of our observations for understanding the physical nature 
of this remarkable system. 

Large Lya nebulae such as SST24 J1434110±331733 are rare. Thus far, most known 
large Lya nebulae are associated with powerful radio galaxies (e.g., McCarthy et al. 1990; 
van Ojik et al. 1997; Pentericci et al. 1997; Dey et al. 1997; Dey 1999; Villar-Martin et al. 
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2002 , 2003; Reuland et al. 2003). Among radio quiet sources, there are two similarly large 
nebulae at z ^ 3.09 (Lya “blobs” LABI and LAB2; Steidel et al. 2000) and one at 2 : ~ 2.38 
(2142—4420 Bl; Francis et al. 1996). These latter systems are associated with large galaxy 
overdensities, and follow-up deeper Lya imaging of these regions has resulted in the discovery 
of numerous other smaller Lya-emitting nebulae (e.g., Palunas et al. 2004; Matsuda et al. 
2004). Similar overdensities have also been found associated with some of the Lya nebulae 
associated with radio sources (e.g.. Keel et al. 1999; Pentericci et al. 2000). 

The rarity of these >100 kpc Lya clouds, their association with powerful AGN and 
galaxy overdensities, and their energetics, all suggest that these regions are the formation 
sites of the most massive galaxies. If so, understanding the physical conditions and energetics 
of these systems can provide important insights into the massive galaxy formation process. 

At present, it is still not clear what powers the nebulae. Is the Lya emission cooling 
radiation of infalling, shocked, pristine gas (Haiman et al. 2000)? Or is it photoionized by 
hot stars or AGN (e.g., Haiman & Rees 2001)? Or is it the cooling of outflowing shocked 
gas (cf. the “superwind model” of Taniguchi & Shioya 2000; Francis et al. 2001)? Much of 
the recent literature on the extended Lya nebulae has focussed on the radio-quiet systems 
(LABI, LAB2 and the sources at z = 2.38), perhaps because the absence of an obvious radio 
source suggested an absence of AGN (and their associated complexity) in these systems. 
However, more recent observations of the LABI system now hnd that it may also contain 
buried AGN (Ghapman et al. 2004). Moreover, there are many similarities between the 
outer, kinematically quieter regions of the Lya halos surrounding the powerful radio sources 
and the radio quiet “blobs”. All this suggests that AGN may exist in all these nebulae: after 
all, if all massive bulges host massive black holes, then perhaps the formation of the active 
nucleus accompanies the formation of the galaxy. 

It is then perhaps not surprising that the known large Lya nebulae reveal such com¬ 
plexity: compact obscured regions, formed galaxies, revealed AGN, all within their extent. 
If these are the largest protogalactic regions, one might expect that both star- and AGN- 
formation could proceed within the cloud. 


4.1. Mass of the Nebula 

The total mass of the nebula can be estimated dynamically (from the velocity prohle of 
the Lya line), under the relatively poor assumptions that the line prohle is not signihcantly 
affected by absorption and scattering. The smooth variation of the central wavelength of 
the line, the relatively constant velocity dispersion, and the overall symmetry of the line 
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all provide some justification to these assumptions. If the velocity prohle observed in the 
inner regions is primarily the result of rotation, then the mass internal to a radius r = 4" 
32 kpc) is M(r < 4") si 0.6 — 1.6 x 10^^ sin“^f Mq, where the large range is due to the 
range in the velocity prohle hts measured from the two different spectroscopic observations. 
The velocity prohle is not well measured at radii larger than ~ 30 kpc. 

Interpreting the line width is also fraught with uncertainty. Absorption by neutral 
hydrogen would result in the measured width being an underestimate of the true width; 
whereas resonant scattering would likely broaden the line. The velocity dispersion of ~ 
400 km s“^ is too large to be due to the thermal motions of the gas. If it is, instead, due to 
the random motions of clouds within the nebula, then we can estimate the dynamical mass 
of the nebula to be M(r < 4") ^ 6 x 10^^ Mq. These dynamical mass estimates derived 
from the line centroid and width are very large. However, they are surprisingly similar to 
the estimates derived for the large Lya halos around high-redshift radio galaxies, where 
dynamical mass estimates range from 0.3 — 10 x 10^^ Mq (Reuland et ah 2003; Villar-Martm 
et ah 2003; van Ojik et ah 1997). We note that in general such dynamical estimates are more 
sensitive to larger masses, since lower mass systems will have gradients that are not resolved 
by typical faint-object spectroscopic observations. 

A separate estimate of the ionized gas mass may be derived by estimating the electron 
density Ue from the net luminosity of the Lya emission line (cf. McCarthy et ah 1990), 
Mion = l.25mpnefvV ^ 1.7 x lO^^/y^(ne/y^/O.O 26 cm“^)M 0 . The hlling factor /y is un¬ 
constrained by the present data, but observations of line-emitting regions in cluster cool¬ 
ing flows suggest hlling factors of order 10“^ to 10“® (e.g., Heckman et ah 1989), which 
would suggest ionized gas masses of Mion ~ 1-7 x lO®(/y/lO“®)M0. If klb = 0.04 and 
f^matter = 0.30, then a comparison with the dynamical mass estimates suggest ionization 
fractions of < 0.013(/y/10“®)°'^(10^^/Mtoi), and a rough upper limit to the hlling factor of 
~ 7.7 X 10“^. The ionized gas fraction is low for a typical hot plasma, suggesting that a 
signihcant fraction of the mass might be contained in colder (neutral), collapsed objects, e.g., 
the galaxies that are forming within the cloud. Although this argument is only speculative 
at present, better measures of the gas density, temperature and hlling factor, perhaps from 
deep rest-frame optical spectroscopy, can help constrain the relative fractions of hot and cold 
gas. 


4.2. Physical State of the Gas and the Nature of the Ionizing Source 

It is difficult to place hrm constraints on the temperature of the gas cloud since at 
present Lya is the only emission line that we detect across the spatial extent of the nebula; 
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He IIA1640, C IVAA1550 and C III]A1909 are only detected in a localized region near the 
center of the nebula, and in addition, C IV and C III] are observed in a red, optically visible 
continuum source which may be associated with the mid-infrared source. There is also weak 
evidence that C IV may be spatially extended between the central region and the MIPS 
source. 

The line ratios of the gas in the central region are not representative of shock ion¬ 
ization: the shock models of Dopita & Sutherland (1996) typically produce much higher 
C IVAA1550/He IIA1640 ratios (e.g., ~ 2 — 20 for shock velocities Ushock ~ 500 — 150 km s“^) 
as compared to the observed value of C IV/He II and tend to overproduce NVA1240 relative 
to what is observed in this source. Higher shock velocities (ushock > 500 km s“^), where 
better matches might be produced, are unlikely given the low velocity widths observed for 
the C IV and He II emission lines. We cannot yet, however, rule out shock ionization in 
other regions of the nebula where Lya is the only detected emission line. 

The He IIA1640 emission line is unusually strong in the central region. Its strength may 
be indicative of low metallicity (and therefore high temperatures) in the ionized gas in this 
region. However, models of star-forming galaxies typically only show He II with rest-frame 
equivalent widths of ^ 15A at very low metallicities (Z < 10“^) and at very young ages 
<2 Myr (Schaerer 2003). Although it is possible we are witnessing this object at such an 
unusual and early epoch, it would be very improbable; it is therefore unlikely that the strong 
He II is due to low metallicity gas. 

Harder ionizing continua, such as those generated by AGN, can also produce strong 
He II; indeed, comparable strengths of C IV and He II are commonly observed in radio 
galaxy spectra and “type 2 quasars” and are generally attributed to photoionization by 
a power-law continuum. The one oddity is the relative weakness of the observed C III] 
emission: most radio galaxy spectra show C III]/C IV ~ 1, whereas we observe a ratio in 
this system of ~ 0.2. The observed ratio can however be produced by a power-law ionizing 
spectrum and is characteristic of a large ionization parameter [/ = Q(H)/47rR^neC ~ 10“^ 
for a /jy oc continuum (Ferland 1981) in this central region. 

A crude constraint on the slope of the ionizing spectrum can be derived from the ratio 
7] of the number of He"*" ionizing photons (~ 3.7 x 10^^ s“^) relative to the number of H 
ionizing photons (r^ 1.6 x 10^® s“^): r/ = (5(He’'')/(5(H) = 0.023 (for a T = 10"^ K gas) (e.g., 
Schaerer 2002). If we relate this ratio to the number of photons at 54 eV relative to that at 
13.6 eV, this ratio suggests a spectral shape of Fjy oc It is intriguing that this slope is 

even steeper than the oc spectrum derived by extrapolation of the mid-infrared SED. 
We have not accounted for Lya absorption in this estimate, which would suggest that the 
ionizing spectrum is even steeper, i.e., a < —2.7. Note that here we have compared the full 
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10" radius aperture measurement of Lya with the spectroscopic measurement of He IIA1640. 
If we compare just the measurements in the same spectroscopic apertures, we derive a ratio 
of f] = 0.21 for T = 10"^ K, which would suggest an extremely hard ionizing spectrum of 
iv . This latter estimate is unlikely to be realistic, since the H+ ionization zone is 

likely to extend over a much larger scale than the zone, but we can take it to be a very 
firm upper limit for the slope of the ionizing continuum. 

If we assume normal HII region electron temperatures of Tg ~ 10^ K in the central 
region, then we can estimate the electron density using the strength of the Lya emission 
in the same aperture used to measure the high ionization lines (see Table 2). The Lya 
luminosity is given as Li^ya = {jhya/jni 3 )nenpfvyhi'i{pa^p ~ 4.05 x where V is 

the volume of the emitting region, Ue and Up are the electron and proton number densities, 
fv is the hlling factor of the ionized gas, and a^^ is the effective recombination coefficient for 
H/5 (Storey & Hummer 1995). In the central region, we find Uefy^ ~ 0.25 cm“^; combined 
with the ionization parameter estimate, this implies a hydrogen ionizing photon flux in the 
central region of Q(H) ^ 3 x fy^'^ s“^ photons s“^. We have no direct measure of 
the electron density, and therefore no hrm constraints on the volume hlling factor of the 
ionized gas. However, if we use the rough upper limit to the hlling factor from the previous 
subsection, then (5(H) > 6 x 10^® photons s“^. 

Using a similar analysis for the total Lya luminosity measured in a 10" radius aperture, 
we hnd ~ 0.026 cm“^ for the average density in the nebula as a whole. For values of 

the hlling factor 10“^ > /y > 10“®, the average electron density in the line-emitting regions 
is in the range 0.8 < Ue < 26 cm“^. 


4.3. Energetics 

The Lya luminosity of the nebula within a 10" radius aperture is ~ 1.67 x 10^"^ erg s“^. 
This is comparable to the luminosities of the other well-known Lya nebulosities, LABI 
(Steidel et ah 2000, 1.4 x 10^"^) and 2142—4420 B1 (Francis et al. 2001; Palunas et al. 2004, 
8 X 10"^^ erg s“^). If it is due to young, hot stars distributed over the face of the nebula in 
a very extended star-forming region, then the luminosity implies a net star-formation rate 
SFRi^ya = LLya/(l-12 X 10^^ erg s“^) ~ 150 Mq yr“^ (Kennicutt 1983). 

It is not clear what source(s) power the Lya emission. No obvious single compact source 
with high UV luminosity is visible within the extent of the nebula. The two galaxies A and 
B do not appear to be plausible candidates for the sources of the photoionizing radiation: 
galaxy B is at the wrong redshift, and galaxy A does not show any strong Lya emission 
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associated with it and its UV luminosity at Arest = ISOOA only corresponds to a modest 
star-formation rate of SFRuy = Lyy/S x 10^'^ erg s“^ Hz“^ 8 Mq yr“^ (e.g., Madau et al. 
1998). 

It is also not clear if the MIPS source is the source of the photoionizing radiation: it is 
faint at rest-frame UV wavelengths and its mid-infrared emission suggests that the source, 
although very luminous, is almost completely obscured from view by dust. It is therefore 
unlikely that short wavelength ionizing photons could escape from the nebula - clearly, 
little or no Lya emission does. Nevertheless, it is intriguing that the estimated far-infrared 
luminosity of the MIPS source is roughly two and a half orders of magnitude larger than the 
Lya luminosity of the entire nebula, and it is possible that some ionizing radiation could 
leak out in directions other than our direct line of sight. We discuss this further below. 

It therefore appears that the Lya emission is powered by sources hidden from direct 
view, perhaps partly by the mid-infrared source, perhaps partly by another obscured or 
mis-directed source lying near the center of the nebula where we also observe the C IV 
and He II emission lines, or by distributed, faint, ionizing sources, or by shocks asso¬ 
ciated with the infall of gas into the region. The dynamical timescale of the nebula is 
given by tdyn ~ |a/R^/2GM ^ 0.37 Gyr. This is much smaller than the age of the Uni¬ 
verse at z = 2.656 (2.4 Gyr in our adopted cosmology), and suggests that we are wit¬ 
nessing the nebula in a young phase. The temperature of the nebular gas is unknown, 
but if we assume a temperature of 10^ K, then the cooling timescale is fcooi ~ ~ 

2.07 X 10^^s(T/10^K)(ne/lcm“^)“^(A(T)/10“^^erg cm^ s“^)“^ (Padmanabhan 2002; Suther¬ 
land & Dopita 1993). For gas between 10"^ — 10^, log(A(T)) can vary from —23.4 to —21.0 
depending on the exact temperature and metallicity, suggesting cooling times which are very 
short (tcooi 1.6 X lO^year), implying that some source of energy (shocks? ionizing photons?) 
is required to maintain the nebula and render it visible. Taken at face value, the large ratio 
of Gyn/Gooi would imply rapid collapse of the cloud into star-forming clumps. 


4.4. The Mid-Infrared Source 

The source SST24 J1434110-I-331733 detected in the Spitzer IRAG and MIPS bands 
lies within the extent of the Lya nebula, roughly 2'.'5 (20 kpc) away from the center. The 
optical spectra show a faint continuum source within the MIPS error circle, just to the 
southeast of the MIPS position. The spectrum of this source appears to be signihcantly 
redder than galaxy A [Fx oc yO-075±o.oo6 MIPS source in the observed optical regime, 

versus Fx oc y-o-076±o.oo4 Galaxy A), but it does show weak associated G III] emission 
and possibly some G IV emission certifying its association with the Lya nebula. The optical 
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spectrum shows no obvious signs of AGN activity (e.g., broad emission lines or strong high- 
ionization lines), but instead shows possible absorption lines similar to those observed in 
young, star-forming Lyman-break galaxies. 

A £t to the SED in the Spitzer bands suggests a power-law spectrum of the form 
Li, ^ 2.7 X 10^^(z//10^^Hz)“^ erg s“^ Hz“^ (figure 11). Rather surprisingly, extrapolating 
this power-law to rest-frame UV wavelengths 1300—ISOOA is consistent with the observed 
spectrum of the faint optical source believed to be associated with the MIPS source, only 
overshooting the observed spectrum by a factor of ~1.3—1.7 (figure 11). This offset is 
perhaps not surprising, given that we are comparing the large aperture total flux densities 
of the Spitzer data with the small aperture (1.5" x 1.0") spectral extraction. 

The mid-infrared source has poorly constrained flux densities at long wavelengths ( > 100 
yum) where much of its luminosity likely emerges. As discussed in section § 3.6, if the SED 
is similar to that of known ULIRGs, the far-infrared luminosity of the source is likely to be 
in the range 1 — 5 x lO^^L©. If this object is a starburst, this luminoisty corresponds to a 
star-formation rate SFR^i ^ 17OO(Lir/1O^^L0)Mq yr“^. If its luminosity-to-mass ratio is 
comparable to ULIRGs (i.e., L/M ^ 100 — 200), the prodigious luminosity implies a total 
mass in dust of about Mdust ~ 0.2 — 5 x IO^Mq if Mdust/Afgas = 0.01. These estimates are 
very coarse, but it is clear that these large bolometric luminosities imply very large masses 
of dust, and that a substantial mass of baryons has already been cycled through stars to 
produce this mass of dust. 

Although the mid-infrared source is extremely energetic and lies within the nebula, it 
is not clear whether it is responsible for photoionizing the nebula. If we assume that the 

power-law represents the intrinsic spectrum of the mid-infrared source, and if it can be 
extrapolated to shorter wavelengths, then we can investigate whether the ionizing photons 
emitted by this source are sufficient to photoionize the nebula and produce the observed 
Lya radiation. Extrapolating the power-law would result in 1.8 x 10®^ ionizing photons 
between 200A and 912A. The luminosity of the Lya nebula requires ~ 10^® ionizing photons. 
The MIPS source can therefore be responsible for at most 18% of the Lya luminosity. This 
is likely a very strong upper limit, since in reality the rest-frame UV spectrum is observed to 
be redder than the extrapolation, and our estimate here does not account for any absorption 
by dust. It is possible, of course, that although the extinction to the source is severe along 
our line of sight, it may still have low optical depth toward the nebula (say, a mis-directed 
quasar). Although plausible, the geometry of the source, taken at face value, does not favor 
such an explanation: the MIPS source is not located at the center of the nebula and instead 
lies ~20 kpc away from the peak of the Lya emission. The symmetry of the Lya emission 
around this peak, the presence of a possible ionizing source near this peak as witnessed by 
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Fig. 11.— The spectral energy distribution of SST24 J1434110+331733, shown here as 
vs V. The rest-frame UV points {v > 10^^ Hz) should be regarded as upper limits, since 
they include continuum contributions from galaxies A and B. The red line at high frequency 
is the spectrum of the optical continuum source closest to the location of the MIPS source 
(shown also in figure 6). The dashed blue line is a £t to the Spitzer data points, oc z/“^. 
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the presence of He II, C IV and C III] emission, and the lack of strong Lya emission cospatial 
with the location of the MIPS source all suggest that it is not primarily responsible for the 
nebula. Further spectroscopy and high spatial resolution imaging of the cloud are necessary 
to carefully account for all the possible ionizing sources within the cloud. 

The power-law SED in the optical and mid-infrared bands exhibited by SST24 J1434110-I-331733 
and its extreme luminosity are more typical of active galaxies and quasars than of star¬ 
forming galaxies (Pier & Krolik 1992; Nenkova et al. 2002, e.g.,). It is therefore likely that 
the infrared luminosity of this source is dominated by an enshrouded AGN, rather than a 
vigorous starburst; starlight from young stars only begins to dominate the flux at rest-frame 
UV wavelengths. In this respect, SST24 J1434110-I-331733 may be very similar to the dust 
enshrouded systems with extreme 24/im-to-optical colors discovered by Houck et al. (2005). 

Spitzer IRS spectroscopy of these extreme systems has established that they are obscured, 
dust-enshrouded AGN at z ^ 2.5 with far-infrared luminosities in the range 6 x 10^^“^^ L© 

(Houck et al. 2005). We speculate that these objects are all extremely massive galaxies 
caught in an early phase of their formation. Since the AGN is currently dominating the 
SED, it is tempting to speculate that the supermassive black hole in these extremely mas¬ 
sive systems form contemporaneously with the assembly of the galaxy, rather than grow 
slowly by subsequent accretion. 


4.5. Galcixy and AGN Formation in the Clond 

In a seminal work, Elmegreen & Lada (1977) proposed that star formation in molecular 
clouds in our Galaxy proceeds through the sequential formation of OB subgroups, with 
each new subgroup driving ionization and shock fronts deeper into the parent molecular 
cloud, stimulating further star formation. Their model was inspired by observations of a few 
Galactic star-forming regions (e.g., W3, M 17 and M 42), where this sequence of events can 
be observed directly: an expanding HII region associated with an OB association located at 
the edge of a molecular cloud. A host of infrared-luminous cores and maser sources lying 
within the molecular cloud and near this edge bear witness to active, ongoing star-formation 
occurring near the location of the shock front, and suggest that the star-formation “front” 
is slowly eating into the cloud. 

It is rather remarkable that the system reported here shows many analogous features: a 
Lyman break star-forming galaxy, which lies just north of a dust-enshrouded active region, 
which lies just north of a large region characterized by Lya cooling radiation. By anal¬ 
ogy to the M 17 region, we postulate that this is a time sequence, with the Lyman break 
galaxy representing the ‘old’, revealed objects, the mid-infrared source representing the site 
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of current star- and AGN-formation, and the central source and Lya nebula representing 
the region where galaxy formation is just beginning. In such an evolutionary scenario, a gas 
cloud, observed by its cooling radiation, would undergo collapse leading to an initial burst 
of star-formation and AGN activity. This immediately enshrouds the active region (through 
the production of dust in the ejecta of first generation supernovae), and then hnally leads to 
a stage when the formed, star-forming galaxy is unveiled. Although highly speculative, this 
scenario would certainly have many merits. It might explain the low stellar masses observed 
in the Lya-bright regions found in the narrow-band surveys (these are the youngest objects 
that have just started forming stars), the very high star-formation rates observed in the 
SGUBA galaxies (these represent the peak star-formation phase), the relatively low star- 
formation rates seen in the Lyman-break galaxies (the early post-starburst phase), and no 
obvious evidence of systems with intermediate metallicities (they are rapidly enshrouded). 

In this case, however, it is very unlikely that the formation events could be triggered 
by an ionization front: over the several kpc scales in question, the only relevant force is 
gravity, and the hypothesis presented here is simply that we may be able to estimate the 
time for the collapse of such a large cloud into its star-forming sub-units by age-dating the 
various sub-units. If, for example, galaxy A is typical of most Lyman-break galaxies at 
2 < z < 3, then we can assume a minimum age of ~ 100 Myr since its last major epoch 
of star-formation (e.g., Papovich et al. 2001). Since its distance from the center of the Lya 
emission is ~ 30 kpc, a naive estimate of the rate at which galaxy formation is proceeding 
through the region is ~ 300 kpc Gyr“^. Better constraints may be derived by detailed age¬ 
dating of the galaxies in the nebula. Our present data preclude such an analysis, but future 
spectroscopic and multiwavelength imaging observations may allow us to date the various 
components more reliably. 

One timescale of relevance here is the formation time of the AGN associated with the 
mid-infrared source. Under the assumption that the bulk of the far-IR luminosity is produced 
by accretion onto a super-massive black hole (SMBH) at the Eddington rate, we estimate 
a central black hole mass of Mbh ~ 2.5 x 10®(L/ij/10^®erg s~^)Mq and an accretion rate of 
M ^ 7.6(L/R/10^®erg s“^)(0.1/?7)Mq yr~^, where rj is the radiative efficiency of the SMBH. 
If the black hole grows at a constant efficiency, then a IO^Mq SMBH can be created by steady 
accretion from a small seed IOMq black hole in about 8 x 10® years, i.e., a seed formation 
redshift of z ^ 4. Our estimates are upper limits for Mbh and the total growth time. For 
more massive seeds or merging of multiple seeds during the formation process, the formation 
time could be signihcantly shorter (cf. Haiman 2004; Yoo & Miralda-Escude 2004; Shapiro 
2005). 



5. Conclusions 


We have discovered a very large Lyo-emitting nebula at a redshift of z = 2.656. This 
is the hrst such object discovered by its associated 24/ini emission. The luminosity of the 
nebula {Li^ya ~ 1-67 x 10"^^ erg s“^ within a 10" radius aperture), its large spatial extent 
(> 100 kpc), the detection of a star-forming “Lyman-break” galaxy and a source of intense 
mid-infrared emission within the nebula all suggest that it is the site of massive galaxy 
formation. In many of these respects, SST24 J1434110-I-331733 is similar to the Lya “blob” 
LABI discovered by Steidel et al. (2000). By analogy to the other large Lya “blobs” known, 
it is probable that this one sits within a large (several Mpc - scale) structure; further deep, 
wide-held Lya imaging of the region should prove fruitful. 

Spectroscopy of the Lya emission reveals that the cloud is kinematically ‘quiet’ at 
least in its inner regions, showing a slowly and monotonically varying velocity centroid with 
position and roughly uniform velocity dispersion across its face. If interpreted as rotation, 
the implied dynamical mass contained within the central 8" is ~ 0.6 — 1.6 x 10^^ Mq. The 
velocity dispersion of the Lya emission line is large (a ^400 km s“^). 

The source of the ionization for the nebula is ambiguous. The geometry of the emitting 
region and the shape of the SED suggest that the mid-infrared source SST24 J1434110-I-331733 
is unlikely to be responsible for more than ~20% of the Lya luminosity, and it is likely that 
the ionization and line radiation are the result of shocks and / or multiple, distributed ion¬ 
izing sources. If the Lya emission is primarily powered by hot, young stars, the total star 
formation rate of the system is ~ 150 Mq yr“^. The detection of localized, weak, narrow 
C IV and He II emission near the center of the nebula suggests that at least portions of the 
nebula may be ionized by a power-law ionizing source or perhaps contain low metallicity gas; 
however, there is clearly 24/im emission from other parts of the system that most likely is due 
to emission from dust, suggesting that the nebula is very inhomogeneous in its properties. 
It is likely that multiple sources of ionization exist within the nebula. 

The mid-infrared source in the nebula is very luminous; if the SED is similar to Arp 220, 
the implied far-infrared luminosity is Lpm ~ 5.8 x 10^^ Lq. The mid-infrared SED is well 
modelled as a power-law oc z/“^, and bears more resemblance to the SEDs of active galaxies 
and quasars than to star-forming galaxies. We suggest that the MIPS source is a forming, 
enshrouded, AGN, perhaps similar to the ~ 2.5 extremely luminous dust-enshrouded AGN 
recently discovered through Spitzer IRS spectroscopy (Houck et al. 2005). Spectroscopy with 
Spitzer IRS will be necessary to determine the relative contributions of a hot dust continuum 
component and PAH emission to the total mid-IR flux, and will help us better constrain the 
relative roles of young stars and AGN in heating the dust. 
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Since the large Lya nebula encompasses a revealed ‘Lyman Break galaxy’ and an 
enhrouded active region, we suggest a scenario in which the two objects represent sites of 
past and present galaxy formation and the Lya emission traces the region of current infall. 
In a manner similar to that observed in nearby HIT regions, the galaxy formation “front” in 
this system is moving gradually to the south (perhaps at a rate of ~ 300 kpc Gyr“^). 

We speculate that Lya nebulae such as this one may be found through their associated 
far-infrared emission. Although this may appear contradictory, it is quite plausible that both 
strong 24/im emisison from dust (powered by enshrouded star-formation and/or accretion ac¬ 
tivity) and strong Lya emission (powered by the cooling of gas in the deep potential well and 
perhaps by photoionizing radiation escaping from the enshrouded regions) are joint signposts 
for the spectacular events associated with the formation of the most massive galaxies. 
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Table 2. Emission-Line Strengths'^ 


Line 

Obs. Wavelength 

(A) 

Aperture 

(") 

Flux 

(erg s“^ cm~^) 

FWHM 

(A) 

Obs. EW 

(A) 

Luminosity'^ 
(erg s-i) 

Lya 

4448 

10.0 radius'^ 

2.89 ±0.15 X 10-1® 

- 

- 

1.67 X loll 

hya 

4444.7 

10.9x1.5 

4.05 ±0.01 X 10-1® 

13.0 

316 

2.33 X lOi® 

Lya 

4445.0 

4.5x1.5 

3.10 ±0.01 X 10-1® 

12.8 

396 

1.79 X lOi® 

NV 

- 

4.5x1.5 

< 5.1 X 10-1® ‘l 

- 

<7.3 


cIV 

5662.0 

4.5x1.5 

4.17 ±0.04 X 10-11 

15.2® 

85 

2.40 X 1012 

Hell 

5997.9 

4.5x1.5 

4.07 ±0.04 X 10-11 

7.3 

60 

2.35 X 1012 

ClII] 

6974.6 

4.5x1.5 

5 ± 1 X 10-1® 

<13 


4.6 X loll 

- 

5081.6 

4.5x1.5 

5.15 ±0.04 X 10-11 

<4 

87 



^All quantities are measured from the PA=19.4° spectrum with the exception of the first entry in the table, 
which is measured from the narrow-band image. 

^ Luminosities are calculated for a cosmology with Hq = 70 km s~^ Mpc~^, = 0.3, Hy = 0.7. 

^Measured from the narrow-band image after subtracting the continuum flux estimated from the G-band 
image 

"^Upper limits quoted are 3cr and assume lines of the same width as the He II feature. 

IV width is measured by modelling the doublet as a single feature. 







